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ABSTRACT 

Context. Substructures - dark lanes and tails - of umbral dots (UDs) were predicted by numerical simulations of magnetoconvection 
and have been detected later in some observations. 

Aims. To provide constraints for realistic theoretical models of sunspot umbrae, we describe the observed properties and evolutionary 
characteristics of UDs (including their substructure) and of other umbral structures. 
f^S, ' Methods. We analyse a 6 h 23 min time series of broadband images of a large umbra in the active region NOAA 10634, acquired 

with the 1-m Swedish Solar Telescope, in the wavelength band around 602 nm. A 43 min part of this series was reconstructed with 
the MFBD method, reaching a spatial resolution of 0714. With the help of image segmentation, feature tracking, and local correlation 
, tracking, we measure brightness, size, lifetime, and horizontal velocities of various umbral structures. 

O !■ Results. Large structures in the umbra - strong and faint light bridges (LBs) and an extended penumbral filament - evolve on time 

scales of hours. Most (90 %) of UDs and bright point-like features in faint LBs split and merge, and their median lifetimes are 3.5 
or 5.7 min, depending on whether the split or merge event is considered as the end of their life. Both UDs and features in faint LBs 
that do not split or merge are clearly smaller (0715) than the average size (0717) of all features. Horizontal motions of umbral bright 
small-scale features are directed either into the umbra or along faint LBs with mean horizontal velocities of 0.34 km s _1 . Features 
faster than 0.4 km s _1 appear mostly at the periphery of the umbra. The motion of peripheral UDs (PUDs) seems to be the continuation 
of the motion of penumbral grains (PCs). The intensity of dark lanes, measured in four bright central UDs (CUDs), is by a factor 0.8 
lower than the peak intensity of CUDs. The width of dark lanes is probably less than the resolution limit 07 14. The characteristic time 
of substructure changes of UDs is ~ 4 min. We observe narrow (0714) bright and dark filaments connected with PUDs. The bright 
filaments are 0.06 / p h brighter than the dark ones. Usually one dark and two bright filaments form a 074 wide tail attached to one 
' PUD, resembling a short dark-cored penumbral filament. 

Conclusions. Our results indicate the similarity between PUDs and PGs located at the tips of bright penumbral filaments. The features 
^-j- . seen in numerical MHD simulations are consistent with our observations of dark lanes in CUDs and tails attached to PUDs. 

£ — . Key words. Sunspots - Sun: photosphere 
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' 1. Introduction Il979h . Another approach, based on the magnetoc onvection in 



. a monolithic thick flux tube, has been proposed by Wei ss et ail 

The study of fine structure in sunspot umbrae is crucial to un- (199Q) Tfl theif model UDs m hot convective plumes that 



derstanding the interactions of plasma with strong magnetic overshoot t0 the sunspot > s photosphere. This approach was ex 

»H ; fields and the physical structure of sunspots. Currently, thanks to tensivel elaborated by means of 3D numerical simulations of 

new ground-based and space instruments image reconstruction magnet oconvection. Recently, Schiissler & Vogler ( 2006) have 

methods, and MHD simulations, our knowledge of fine-structure , ,,, . . ■ .. , .. c e 

. . . , , ' £, . shown that the energy transport in the vertical magnetic field of 

elements is showing considerable progress This concerns urn- ^ umbra ^ dominated b nonstationary nalT0W plumes f ris . 

bral dots (UDs , the small isolated bright features observed in ■ hot lasma wkh adjacent downflows . The magne tic field is 

the umbra, and light bridges (LBs), the elongated bright features . , ■ , . , e . u , .... • -,, 

' , » , , , , , strongly reduced in the upper layers of the plumes. At the visible 

crossing the umbra. At the umbra-penumbra boundary penum- surface (r = 1} ^ lumefj m geen ag brfght features whh a 

bral grains (PGs) located at the tips of bright penumbral ffla- ^ size of 20Q _ 30Q km and nfetime of the order of 3Q ^ 

ments sometimes protrude deep into the umbra. Usually, UDs , ,■ TTT> . e ■ , . , TTrv , , . , 

F " , , , , , resembling UDs. Most of the simulated UDs have an elongated 

are divided into two groups peripheral umbral dots (PUDs) and sh ^ a central ^ ^ some j UDs show a threefold 

central um bral dots (CUDs), according to their position inside ^ kne The ^ lanes ^ a conse e of enhanced density 

* e umh ™ d Grossmann-Doerthetal.| | l986|) . For definitions and and d at ^ of ^ j SQ ^ ^ continuum forma . 

basic characteristic s of the sunspot fir^sttuctures, we refer to tion , . shifted dg tQ coolej . k of ^ atmo here . 

the monograph by |Thomas & Weiss| 42008|). iRempel et alJ J2009h extended the simulation, introducing a 2D 

Two basic approaches have been p ror^ed to explain the na- (&m m et ic field fanning out with height into a 3D convec- 

tU /t UDs. According to| Choudhuri| (lH86j), UDs are columns tiye media Th obtained a sun t . like structure surrounded by 

of hot field-free gas intruding into a cluster of thin magnetic photospheric granu i ation . i n the umbra of the "sunspot", UDs 

flux tubes that form the m agnetic structure of a sunspot ([Parkeij similar tQ thos£ simulated by ISchiissler & Voglel ([M were 
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netic field (with respect to the vertical), UDs appeared with two 
attached bright narrow tails pointing toward the penumbra. Even 
farther from the centre of the umbra structures developed re- 
sembling PGs followed by da rk-cored penumbral filaments (see 
Fig. 1 of iRempel et alj|2009t) . 

The most recent ext ensive imaging observa tions of UDs 
have been published by Riethm uller et al.l d200 8b). They ob- 
served thousands of UDs in a 106 min long time series of im- 
ages acquired in the TiO band head (705.7 nm) with the 1-m 
Swedish Solar Telescope (SST). Applying a feature tracking al- 
gorithm, they measured sizes, brightnesses, lifetimes, trajecto- 
ries, and horizontal velocities. The majority of UDs was spa- 
tially resolved with typical diameter of 0731. A difference in 
the evolution of brightness is found between CUDs and PUDs, 
pointing to a different origin of these types of UDs. The be- 
ha viour of CUDs shows bette r agreement with the simulations 
of lSchiissler& Voglerl (120061) than that of PUDs. A new sub- 
set of PUDs, mobile UDs with travel distances over 1", was in- 
troduced. Mobile UDs are born close to the umbra-penumbra 
boundary, they are brighter, live longer, and move faster (680 
m s -1 ) than the rest of UDs (410 m s" 1 ). Pro perties of UDs 



were also recently studied by iKitai et all d2007l) in three series 
of broad-band images acquired by Hinode/SOT. 

Recent spectroscopic and spectropolarimetric observations 
have provided new information about line-of-sight (LOS) veloc- 
ities and magnetic fields. Analysing Dopplergr ams obtained at 
the Dunn Solar Telescope, iBharti etaU d2007al) found upflows 
of 400 m s surrounded by downfiows of 300 m s _1 in large 
(075) UDs. The discovery of downfiows strongly supports the 
magnetoconvective origin of UDs. Full Stokes spectra of 21 
C UDs and 30 PUDs , observ ed by Hinode/SOT, were inverted 
by iRiethmulle r et all d2008a). At the continuum level (r = 1), 
they find a magnetic field weakening of 480 G in CUDs and 510 
G in PUDs. No significant LOS velocities have been detected 
in CUDs, but PUDs are connected with upflows of 800 m s _1 . 
The physical differences between CUD s and PUDs have been 
confirmed by Sobotka & Jurcak ( 2009), using a 3 h long time 
series of full-Stokes spectral scans obtained with Hinode/SOT. 
In the range of optical depths 0.3 < r < 0.6, CUDs do not show 
any excess of LOS velocity and magnetic field inclination with 
respect to the surrounding umbra, while upflows of 400 m s _1 
and a more horizontal magnetic field are observed in PUDs. It 
seems that PUDs are more like PGs than CUDs. Both PUDs 
and PGs appear in regions with a weaker and more horizontal 
magnetic field and their formation height reaches the low pho- 
tosphere, while CUDs appear in regions with stronger and more 
vertical magnetic fields, and they are formed too deep to detect 
upflows and changes in magnetic field inclination. 

Dark lanes in UDs, predicted in the sim ulations by 
ISchussler & VogleJ d2006l) . have been observed by IBharti et al.l 
d2007bl) in large (>1") bright UD-like features seen in 
Hinode/SOT G-band images with diffraction-limited spatial res- 
olution of 0722. After inverting a nearly contempora neous 
Hinode/SOT spectropolarimetric scan, Bha rti et al.l {2009) con- 
clude that the dark lanes correspond to a temperature deficit lo- 
cated above t = 0.1, which, together with observed upflows and 
magnetic field reduction in UDs, supports the results of theoreti- 
cal simulat ions. In G - band i mages acquired with the Dunn Solar 
Telescope, Rimme lei (120081) observed dark dots or lanes within 
bright UDs. Their size is typically close to the diffraction limit of 
07 12. Close to the umbral border, Rimmele finds elongated dark 
structures associated with bright PUDs moving into the umbra. 
These dark structures have the appearance of a short dark tail 



pointing toward the penumbra, and they resemble dark-cored fil- 
aments originating at PGs. 

In this paper we describe the fine structure of a large and 
dark umbra, observed for more than 6 hours with the spatial 
resolution 0714 in a large number of broad-band images. First 
we describe the observed long-term characteristics and evolu- 
tion of the umbra, LBs, and long penumbral filaments extending 
into the umbra. Then we centre on some statistical properties of 
umbral features. Some preliminar y results were published pre- 
viously in two conf erence papers (So botka & Pu schmann 2007; 
Sobo tka et al.l|2008l hereafter Papers I and II). Finally we focus 
on the substructures of UDs - dark lanes in CUDs and fine bright 
and dark filaments composing tails associated with PUDs. 

2. Observations and data processing 

The large leading sunspot in NOAA 10634 was observed on 18 
June 2004 from 07:43 to 15:30 UT with the 1-m SST, La Palma 
dScharmer et al.ll2003l) . equipped with adaptive optics. The spot, 
located near the disk centre (N13, E12), was growing to the max- 
imum phase of evolution. The broad-band images were acquired 
simultaneously in three wavelength bands: blue (450.8 ±0.5) nm, 
red (602.0 ± 1.3) nm, and G-band (430.89 ± 0.6) nm. Exposure 
times were 11-14 ms and the pixel size was 070405 x 070405. 
The images were acquired in the frame-selection mode during 
most of the observing period. The selection interval was 20 s 
and the average time difference between two frames (including 
the time necessary to store the images) is 20.1 s. From 08:40 to 
09:30 UT, all frames (cadence 4-5 s _1 ) were stored to be used 
later for image reconstruction by MFBD (the special treatment 
of these data is described in this section later on). A post-facto 
frame selection was applied to these frames to obtain a homoge- 
neous time series during the full observing period with a tempo- 
ral resolution of 20.1 s. 

After dark- and flatfield corrections, the intensities in each 
image were normalised with respect to the average intensity 
7 p h of a quiet granulation area near the sunspot to compensate 
for changes of transparency and/or exposure time. Photospheric 
and penumbral stray light, which contaminates the umbra due 
to the scattering in the terrestrial atmosphere and in the in- 
strume nt, was eliminated using the method by Martinez Pillet 
(119921) . Parameters of the scattering were determined for each 
wavelength from the shapes of photometric profiles across the 
solar limb. Images of the limb were recorded immediately after 
the observation. The level of stray light, expected to originate 
mostly in the instrument, was 8.5 % in blue, 6.5 % in red, and 
8.7 % in G-band. A deconvolution of the instrumental profile of 
the diffraction-limited 1-m telescope and a simultaneous noise 
filtering was applied to the frames of the series, using Wiener 
filters with noise suppression starting at O/ 1 1 (blue), 0714 (red), 
and 0713 (G-band). Regarding the correction of wavefront aber- 
rations done by the adaptive optics, these values characterise the 
spatial resolution in the best frames. Then, the image rotation 
was compensated and the frames were aligned and de-stretched. 
Finally, a subsonic filter with a cutoff at 4 km s _1 was applied 
to remove a residual jitter due to the seeing. The beginnings and 
ends of the series were removed because of insufficient image 
quality and apodisation effect. The resulting three time series 
in blue, red, and G-band wavelength ranges span from 08:05 to 
14:28 UT, i.e., 6 h 23 min. The G-band series is not suitable 
for the study of umbral structures, because the darkest parts of 
the umbra are underexposed. The red continuum series with the 
best signal-to-noise ratio (S /N = 200) was selected for further 
analysis. 
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As mentioned above, a part of the time series (from 08:40 
to 09:30 UT) has been treated with special care. During this 
50 min time period the seeing was not extremely good but 
very stable, offering a good opportunity for an efficient image 
restoration. We stored all the images (a total of 12000 frames 
in each blue and red band) acquired at a frame rate of 4-5 im- 
ages s _1 . We applied the post-facto multi-frame blind deconvo- 
lution (MFBD) image restora tion algorithm developed and im- 
plemented by Lofdahl (2002]) on the recorded dark- and fiat- 
field corrected images. This reconstruction method uses multi- 
ple frames to retrieve the residual aberrations affecting the im- 
ages after the partial correction performed by the adaptive op- 
tics and thus further improves the image quality. In the present 
case we combined sets of 50 sequential images to obtain sin- 
gle restored ones. Thus, the resulting sequences of the red and 
blue bands last 50 min and consist of 236 reconstructed images 
with a temporal cadence of 12.6 s each. This sampling interval 
has been chosen such that the amount of images grouped for 
the restoration is as large as possible without losing important 
physical information about the temporal evolution of our sci- 
ence target. Moreover, the reconstruction is performed in small 
portions of the images where one can consider a common char- 
acteristic aberration (isoplanatic patches). Overlapping patches 
of 128x128 pix (5" 184x5" 184) have been extracted such that 
at the end of their restoration the whole field of view (FOV) 
can be recovered by mosaicking them without any discontinu- 
ity. Proceeding in this way, the homogeneity and image quality 
along the time series and in the whole FOV can be improved 
up to close to the diffraction limit of the telescope. The result- 
ing reconstructed time series have been treated with the standard 
procedures, i.e., image de-rotation, image alignment, and sub- 
sonic filtering (again with a cutoff of 4 km s _1 ). Some frames 
of these series were lost due to apodisation, so that the resulting 
43 min long blue- and red-band series contain 204 frames each 
and span from 08:43 to 09:26 UT. 

The horizontal motions of umbral features were stud- 
ied using the methods o f local correlation tracking (LCT, 
iNovember & S imon 1988) and feature tracking. LCT provides 
a time-averaged horizontal velocity field of all structures in the 
FOV. The FWHM of the Gaussian tracking window was set 
to 0732 and the temporal integration was made over 1 h and 
2 h. There was no substantial difference between the veloc- 
ity fields obtained with these two integration times. The fea- 
ture tracking provides time records of positions, intensities, and 
sizes for each individual featu re. We used the algorithm devel- 
oped by Sobot ka et alj d 19971) that has recently been improved 
bv lHamedivafal d2008l) . This algorithm makes it possible to dis- 
tinguish between features that split into two or more parts, merge 
with another feature, or preserve their integrity during their life. 

3. Results 

3.1. Long-term evolution of umbral structures 

The evolution of the umbra was inspected visually in the movie 
composed of the red-band images taken from 08:05 to 14:28 
UT. One of the images, taken at 12:21 UT, is shown in Fig.Q] 
The umbra was very stable during our observation. Its area 
A = 345 arcsec 2 , corresponding to the effective diameter d e g = 
yfAAjn = 21", was constant within this period. The area was 
defined by the intensity threshold 0.5 / p h. The observed umbra 
is very dark. After stray-light correction, the minimum inten- 
sity in the umbra is 0.09 7 p h (0.05 7 p h in the blue band), re- 
maining constant with time. The whole umbra shows a slowly 




Fig. 1. Red-band image of the umbra taken at 12:21 UT. 1 - 
strong light bridge; 2,3 - faint light bridges; 4 - long penumbral 
filament. 



evolving grainy and/or filamentary low-contrast "background" 
pattern with rms fluctuations of only 0.01-0.02 7 p h (Paper I). 
Bright inward-moving UDs travel in this background pattern. 
Even a dark nucleus, the less intense part of the umbra, where the 
strongest magnetic field is expected, shows a very faint grainy 
and filamentary structure. Identical patterns were detected in the 
red and blue bands, proving that we observe real structures, so 
the umbral background does not have the diffuse character ex- 
pected in the past. 

A strong LB, marked as "1" in Fig.Q] divided the umbra into 
two umbral cores. One part of the bridge (top right in Fig. Q]) 
was composed of extended penumbral filaments with dark cores, 
while the other part had a granular structure with a central dark 
lane. The dark lane was not connected to any of the dark cores 
in the penumbral filaments. After 3 h 5 min of the observation 
(at 11:10 UT), the granular part split into several fragments that 
weakened gradually, until the minimum brightness was reached 
at 13:35 UT (see snapshots of evolution in Fig.|2j top row). Then 
the granular part started to recover. Bright grains moved along 
the filamentary part of the LB toward the granular one with typ- 
ical velocities of 0.8-1 km s _I (LCT) all the time. 

The long penumbral filament "4", extending deep into the 
umbra, passed several evolutionary phases. It had been present 
since the beginning of the observation. Starting from 08:22 UT, 
its length increased with an average speed of 0.6 km s _1 dur- 
ing 40 min. After that growing phase, it remained stable for 50 
min and then it started to grow again. The second growing phase 
took approximately 30 min. In the next 20 min the filament sep- 
arated from the penumbra and then, after a 30 min long decay 
phase, it disappeared at 11:10 UT. The bright head of the fila- 
ment converted into three PUDs that moved toward LB 2. Three 
snapshots displayed in the bottom row of Fig. [2] show this fila- 
ment in the second growing phase, during its separation from the 
penumbra, and in the decay phase. At the same time (11:10 UT), 
a new long penumbral filament developed at the distance of 173 



4 



M. Sobotka and K. G. Puschmann: Morphology and evolution of umbral dots 



Aid 

£ 




f L 

W v . « 

IT 13:35 


10:40 


i 


*■ 

10:58 



Fig. 2. Snapshots of evolution of the strong light bridge marked 
as "1" in Fig.Q](top row) and of the long penumbral filament "4" 
(bottom row). 

from the original one. It showed the same phases of evolution 
at very similar time scales - growth (50 min, speed 0.6 km s _1 ), 
stability (50 min), second growth (30 min), separation (20 min), 
and decay (30 min). It disappeared at 14: 1 1 UT and its total life- 
time was 3 h. During the existence of both filaments, bright PGs 
moved along them with speeds up to 1 .2 km s _1 (LCT) and many 
PUDs, moving toward LB 2, detached from the heads of the fil- 
aments. 

Faint granular LBs "2" and "3" (Fig. [U started to form at 
10:05 and 09:15 UT, respectively. They developed into well- 
distinguishable structures with central dark lanes around 1 1 :45 
UT. During the formation period, the motion (0.2-0.3 km s _1 , 
LCT) of bright features along LB 3 and along the upper part of 
LB 2 was oriented into the umbra. In the lower part of LB 2, 
an opposite motion toward the penumbra (0.3 km s _1 ) was seen. 
When the formation was finished, the opposite motion of 0.4 
km s" 1 prevailed in LB 2 (see Paper I). Such opposite motion s 
along two adjacent LBs were also observed by iRimmelel (12008). 
This LB started to decay at 13:55 UT, while LB 3 remained 
stable till the end of the observation. The opposite motion of 
bright features in LB 2 might be related to a horizontal velocity 
field connected with two long penumbral filaments that devel- 
oped subsequently at the position "4" in Fig. Q] An example of 
an LCT map of horizontal velocities in the umbra is shown in 
Fig. [3] This map is calculated for the period from 12:15 to 14:12 
UT. 

Bright inward-moving PUDs are often followed by faint 
bright filaments, resembling penumbral filaments attached to 
PGs, and/or by dark "wakes", resembling dark cores of penum- 
bral filament s (Paper I) . Simi lar structures were reported inde- 
pendently by IRimmelel (12008b . The LCT flow map (see Fig. [3) 
shows that the horizontal motion of PUDs (0.5 km s ) seems to 
be the continuation of the inward-directed motion field of PGs 
in the inner penumbra. 

3.2. Statistical properties of umbral features 

In the following we focus on some detailed statistical results ob- 
tained for the best part of the unreconstructed series in the red 
band, taken in the frame-selection mode from 12:15 to 14:12 
UT. Movies covering 90 minutes of the time series from 12:28 
to 13:58 UT can be found at http : //www . asu . c as . cz/~ sdsa/ 
gallery-astropictures . html. 

The feature-tracking technique was applied to this 117 min 
long unreconstructed time series. The bright features were sepa- 
rated from the background using the low-noise curvature deter- 




Fig. 3. LCT map of horizontal velocities, calculated with the 
tracking window of 0732 and averaged over a period of 1 17 min- 
utes, the best part of the unreconstructed series in the red band, 
taken in the frame-selection mode from 12:15 to 14:12 UT. The 
length of the black horizontal bar at the bottom left corner rep- 
resents 1 km s . 

mination algorithm (Hamedivafa 2008). Only features (UDs and 
dot-like features in faint LBs) brighter by 0.04 7 p h than the sur- 
rounding umbra, larger than C'14 (resolution limit), and living 
longer than 2 min were taken into account. The resulting sample 
contains 3099 UDs and 1357 dot-like features in LBs. 

The majority (90 %) of observed features split or merge 
with other features (see Paper II), thus it is necessary to take 
into account these effects. One possibility is to finish the time 
record of a feature when it splits or merges. In this case, the 
mean (median) lifetimes of UDs are 4.5 min (3.7 min) and the 
mean (median) lifetimes of the features in faint LBs are 4.4 min 
(3.4 min). On the other hand, split and merge events caused by 
seeing effects and random coincidences may spuriously shorten 
the lifetime of many UDs. The used feature tracking algorithm 
makes it possible to connect the records of features before and 
after the split/merge event according to the following rules: (a) 
if a feature splits, its successor is the brightest fragment; (b) if 
two or more features merge, the resulting feature is connected 
to its longest-lived predecessor. In this case, the mean (median) 
lifetimes of UDs are 9.1 min (5.7 min) and the mean (median) 
lifetimes of the features in faint LBs are 9.8 min (5.7 min). 
Independently of the applied method, the histograms of life- 
times show a monotonous decrease from shortes t to lo ngest 
lifetimes as reported pre viously by ISobotka et al.l d 1997b and 
Riethmuller et al. I (l2008bl) . 

Both UDs and features in faint LBs 2 and 3 that do not 
split or merge are smaller than the average size of all features 
(Paper II). Histograms of effective diameters calculated from 
time-averaged areas of all 3099 UDs and 1357 bright features 
in LBs and of 337 UDs and 1 14 bright features in LBs that did 
not split or merge are plotted in Fig. [4] top row. It can be seen that 
UDs and bright features in LBs do not differ in size. However, 
the non-split/merge features have narrower range of sizes and 
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Fig. 4. Histograms of time-averaged effective diameters (top) 
and time-averaged horizontal velocity magnitudes (bottom), de- 
rived from feature tracking. Dotted line - non-split/merge fea- 
tures, solid line - all features. 



smaller mean size (0715) than all features in the sample (with 
the mean size of 07 17). Larger diameters of splitting or merging 
features suggest that these features may be composed to a large 
extent of small unresolved objects. 

The feature tracking was also used to measure the horizon- 
tal velocities of small-scale umbral features. Motions of UDs 
are mostly directed into the umbra, while the features in faint 
LBs 2 and 3 move predominantly along the bridges. Histograms 
of time-averaged velocity magnitudes are plotted in Fig. [4] bot- 
tom row. We see that the velocity magnitudes are similar for UDs 
and features in LBs and do not depend on split or merge of the 
features. All types of features show a mean velocity around 0.34 
km s -1 . A map of horizontal velocity vectors (not shown here) 
indicates that features faster than 0.4 km s _I appear mostly at 
the periphery of the umbra and in faint LBs, where the magnetic 
field is expected to be weaker and more horizontal (see Paper II). 

The time-averaged values of brightness, size, and horizontal 
velocity of small-scale umbral features are completely uncorre- 
lated. However, during the evolution of an individual feature, 
there is a correlation between the brightness excess AI with re- 
spect to the surrounding umbra and the area of the feature (Paper 
II). An example of the evolution of a long-lived CUD is shown 
in the left panel of Fig. [5] The correlation between AI and area 
was 0.68 in this case. Coefficients of correlation between AI and 
area evolution were calculated for the non-split/merge features 
(337 UDs and 1 14 bright features in faint LBs 2 and 3) and the 
histograms are plotted in the right panel of Fig. [5] The temporal 
brightness and size variations are positively correlated in most 
cases - 50 % of features have correlation coefficients higher 
than 0.48 and the typical value of a correlation is 0.60. The fea- 
tures first grow and brighten and then get s maller and fainter, in 
agree ment with theoretical predictions (e.g. Schussle r & Voglerl 
2006). For the other features that split or merge, the A/-to-area 
correlation is much weaker but still detectable. The histograms 
(not shown here) are asymmetric with prevailing positive values 
and have flat maxima around 0.40. For 50 % of split/merge fea- 
tures, the correlation coefficients are higher than 0.22. The rea- 
son for the weak correlation is that the split/merge events intro- 
duce substantial uncorrected changes in brightnesses and sizes 
of the features. 
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Fig. 5. Left - example of area (solid) and brightness (dashed) 
evolution of a central umbral dot. Right - histograms of the cor- 
relation between brightness and area during the evolution of in- 
dividual non-split/merge features. Solid line - umbral dots, dot- 
ted line - bright features in faint light bridges. 




Fig. 6. MFBD reconstructed red-band image of the umbra taken 
at 08:59 UT. This image is enhanced by the unsharp masking 
technique for a better visualisation. The upper rectangle contains 
a group of central umbral dots with dark lanes, the lower one a 
field of peripheral umbral dots with tails. 

3.3. Dark lanes in central umbral dots 

A group of bright (0.3-0.5 7 ph ), large (074), and long-lived CUDs 
was present in the central part of the umbra from 08:30 to 10:05 
UT (see Fig. [6] upper box). This group moved slowly in the di- 
rection of LB 2. Dark lan es, predicted by theoretical simulations 
(ISchussler & Voglerl2 006). are observed in several CUDs of this 
group. Obviously, because of the limited spatial resolution, we 
are able to observe dark lanes only in bright and large UDs. 
According to simulations, such UDs should show multifold dark 
lanes, separating them into several parts. In real observations, it 
is difficult to distinguish such UDs from conglomerates of unre- 
lated small UDs. However, the integrity of the observed CUDs 
during the time period of about one hour indicates that we indeed 
see CUDs with dark lanes. 

The 43 min series of MFBD reconstructed red-band frames 
from 08:43 to 09:26 UT with a temporal resolution 12.6 s is 
used to study the substructures in the four brightest CUDs in de- 
tail. Equivalent substructures are also seen in the reconstructed 
blue-band frames, but with lower signal-to-noise ratio. Three- 
or fourfold dark lanes separate CUDs into 3-4 parts. The sub- 
structures vary with time (the typical time scale is about 3 min) 
and the dark lanes disappear and reappear during the evolution 
of CUDs. These changes may be caused by residual seeing vari- 
ations. 

The measurement of brightness, size, and temporal changes 
of four CUDs, their parts, and dark lanes was made in a se- 
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Fig. 7. Central umbral dots with dark lanes. The field corre- 
sponds to the upper box in Fig. [6] Left - original image, right 
- corresponding segmented image. Umbral dots under study are 
marked by numbers. 



Table 1. Observed parameters of central umbral dots with dark 
lanes (see text for the explanation of symbols). 



CUD 


Aid 


h 


/di 


d 


t [min] 


f d i [min] 


1 


0.54 


0.14 


0.43 


0741 


80 


35 


2 


0.36 


0.15 


0.28 


0740 


75 


10 


3 


0.43 


0.18 


0.32 


tt' 42 


70 


14 


4 


0.43 


0.16 


0.37 


0.' 44 


50 


14 



ries of segmented images of the field delimited by the upper 
box in Fig. [6] The segmentation was made in two steps. First, 
CUDs were se parated from the backgr ound using the multilevel 
segmentation (Bovelet & Wiehrll200ll T. Roudier 2003, private 
communication). Then the dark lanes and parts of CUDs were 
identifie d by means of the l ow-noise curvature determination al- 
gorithm dHamedivafal2008b . An example of the original and seg- 
mented image is displayed in Fig. [7] 

The basic observed parameters of CUDs \-A are summarised 
in Table Q] The table contains the values of CUD peak intensity 
7 u d, intensity 4 of the surrounding umbral background, intensity 
/di of the dark lanes, and diameter d. All the intensities are in 
units of /ph. The values of intensities and diameters are the aver- 
ages of measurements in four frames for each CUD showing the 
best-resolved substructure. The lifetime t of CUDs is measured 
in the 6 h 23 min long time series and t& is the total time for 
which the dark lanes were observed during the 43 min long series 
of reconstructed frames. The peak-to-background ratio I u d/h f° r 
the brightest CUD 1 is 3.9, for CUDs 2-4 is approximately 2.5. 
The average ratio is 0.8. Since the dark lanes are very 

narrow, this value has to be considered as an upper limit. 

The feature-tracking algorithm was applied to the 43 min se- 
ries of segmented images to measure the lifetimes and sizes of 
the individual parts of CUDs, separated by the dark lanes. The 
mean and median lifetimes are 9 min and 4 min, respectively. 
The median value is comparable to the characteristic time of 
temporal variations of the substructures (3 min) estimated by the 
visual inspection of the movie. Residual seeing fluctuations may 
be the cause of these changes. Both the mean and median sizes 
of the parts of CUDs are 0717. The apparent width of the dark 
lanes measured in the segmented images is 0708-0712 (2-3 pix- 
els), i.e., smaller than our resolution limit 0714. Nevertheless, 
because the dark lanes are easily detected in the restored images 
(see Fig. [7] left panel), we may expect that they are not much 
narrower than 07 14. 




Fig. 8. Narrow bright and dark filaments (tails) connected to pe- 
ripheral umbral dots. Numbers label examples of: 1 - bright fil- 
aments, 2 - dark filaments, 3 - tails composed of two bright and 
one dark filament. The image is enhanced by the unsharp mask- 
ing technique for a better visualisation. The field corresponds to 
the lower box in Fig. [6] 



3.4. Tails of peripheral umbral dots 

In Paper I we mentioned that narrow bright or dark tails are 
often attached to PUDs. In the 43 min series of MFBD recon- 
structed red- and blue-band images, a well-resolved pattern of 
radially oriented narrow bright and dark filaments is observed 
in the peripheral parts of the umbra (see Fig. [6]). These fila- 
ments are not connected to the penumbra. The rms contrast of 
the pattern is only 0.02-0.03 7 p h in the red-band images. Most 
of the filaments are attached to PUDs, giving the appearance of 
come t-like tails poi nting in the direction toward the penumbra 
(cf. lRimmelel l2008T). The presence of the attached filaments is 
independent of the brightness of PUDs. An example is shown 
in Fig. [8] a magnification of the lower box in Fig. [6] The image 
is enhanced by unsharp masking for a better visualisation of the 
small-scale structure. Examples of bright and dark filaments are 
labelled by numbers 1 and 2, respectively. 

The basic observed parameters of the bright and dark fila- 
ments were measured in the best images of the reconstructed 
red-band series. The typical width of the filaments is at or be- 
low the resolution limit 0714. Their length ranges from 074 to 
1". The intensity of the bright filaments is approximately 0.06 
7 p h higher than the intensity of the dark ones. The filaments live 
5-20 min, what is comparable with the lifetime of UDs. 

In many cases, two bright filaments with a dark one between 
them are attached to one PUD, forming a composite tail approx- 
imately 074 wide, similar to a short dark-cored penumbral fil- 
ament. Such dark-cored tails are seen in Fig. [8] and also in the 
right part of Fig. [6] at the periphery of the umbra and near the 
strong light bridge. Sometimes we observe that one bright fila- 
ment is missing, probably because it is too faint and narrow to 
be detected, sometimes the dark-cored tails overlap. For this rea- 
son we expect that the dark-cored tails of PUDs are in fact more 
numerous than we observe. Such composite tails also appear as 
a result of the theoretical simulations of magnetoconvection in 
an inclined magneti c field. They are well visible in Fig. 1 of 
Remp eTet al.l (|2()09) at the periphery of the simulated umbra. 
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4. Discussion and conclusions 

We analysed the 6 h 23 min time series of broadband im- 
ages (A = 602.0 ±1.3 nm) of a large (d e s = 21") and dark 
(7min = 0.09 7 p h) umbra. A 43 min part of this series wa s re- 
constructed by means of the MFBD method (Lofdahl 2003) and 
used to study the substructures of UDs - dark lanes in CUDs and 
narrow tails attached to PUDs. 

Two faint LBs with dark central channels (LB 2 and LB 3, 
see Fig. Q]) were formed in the umbra. The formation phase 
took about 2 h, followed by the phase of stability. Small (0717) 
bright dot-like features moved along these LBs. The dark cen- 
tral channels in faint LBs may have a similar physical origin like 
the dark lanes in UDs, i.e., the enhanced density at the top of 
ascending hot gas dRimmeldl2008t iGiordano et alj[2008). Two 
bright penumbral filaments extending deep into the umbra devel- 
oped subsequently at nearby positions ("4" in Fig.[T|). Each one 
showed two phases of growth (50 and 30 min), then a separation 
from the penumbra (20 min), and a final decay into several PUDs 
(30 min). The total lifetime of one penumbral filament was 3 h. 

Lifetimes, sizes, and horizontal velocities of UDs and dot- 
like features in faint LBs were measured using a feature track- 
ing technique applied to the 117 min long best part of the un- 
reconstructed time series. The majority (90 %) of all observed 
features split or merge with other features. The median lifetime 
of both types of features is approximately 3.5 min if the life of 
a feature is expected to finish when it splits or merges. If the 
split/merge events are not considered to be the end of life (ac- 
cording to the rules described in Section 13.21 ), the median life- 
time is 5.7 min. Both UDs and features in faint LBs that do not 
split or merge are clearly smaller (0715) than the average size 
(0717) of all features. This means that features showing the split 
or merge events may be composed to a large extent of small 
unresolved objects. The obtained sizes are much s maller com- 
pared to 073 1 reported bv lRiethmuller et alj d2008bl) . Horizontal 
motions of umbral bright small-scale features are directed either 
into the umbra or along faint LBs. Time-averaged velocity mag- 
nitudes (0.34 km s~') are similar for UDs and features in faint 
LBs and do not depend on the split or merging of the features. 
Features faster than 0.4 km s _1 mostly appear at the periphery of 
the umbra and in faint LBs, where the magnetic field is expected 
to be weaker and more horizontal. Although the time-averaged 
values of brightness, size, and horizontal velocity of features are 
uncorrelated, during the evolution of individual non-split/merge 
features the temporal brightness variations are positively corre- 
lated with the size variations in most cases. 

For the first time, we tried to measure the detailed photo- 
metric parameters of CUDs with dark lanes in the series of re- 
constructed images. For this purpose we selected a group of 
four bright, large, and long-lived CUDs in the centre of the um- 
bra. They were brighter by a factor of 2.5 (3.9 in one case) 
than the surrounding background. The three- or fourfold dark 
lanes, observed with the resolution of 0714, correspond to nar- 
row - below the resolution limit - 20 % deep depressions in 
the intensity profiles of CUDs. This result has to be taken with 
care, because of the strong contamination by stray light from 
the bright parts of CUDs. The substructures vary strongly with 
time. The median lifetime of the bright parts of CUDs is 4 
min, and the dark lanes disappear and reappear. It is a ques- 
tion of whether these cha nges, not seen in the simulations by 
Schiissler & Vogler ( 2006), might be caused by the residual see- 
ing variations in the reconstructed images. The average size and 
lifetime of the four selected CUDs with dark lanes (0742 and 70 
min) are comparable to the corresponding parameters of UDs ap- 



pearing in the simulations (073-074 and 30 min). However, such 
large and long-lived CUDs are not frequent in umbrae. The ma- 
jority of UDs are sm aller (072-073) and live m uch shorter time 
(3-9 min) - see also lRiethmuller et al.l d2008bl) . The future the- 
oretical simulations should reproduce also such small and short- 
lived UDs. 

We confirm the bright and dark filamentary tails connected 
to PU Ds and directed to ward the penumbra, reported in Paper I 
and bv lRimmeld d2008l) . In fact, such substructures are attached 
to nearly all PUDs, independent of their brightness. We newly 
observe that many tails of PUDs are composed of two narrow 
bright and one dark parallel filaments, forming a structure 074 
wide that resembles a short dark-cored penumbral filament. The 
length of this structure is 074— 1". The lifetime of such dark- 
cored tails is comparable to the lifetime of UDs. It can be ex- 
pected that the origin of central dark filaments in the tails is anal- 
ogous to the origin o f the dark lanes in UDs a nd dark cores in 
penumbral filaments. Sobotka & Jurcak] d2009l) found that, con- 
cerning the physical characteristics, PUDs are similar to PGs lo- 
cated at the tips of bright penumbral filaments. The discovery 
of the dark-cored tails attached to PUDs strongly supports this 
similarity. Such structures also app ear as a result o f the r adiative 
MHD simulation of a sunspot by iRempel et alj d2009t) . These 
authors claim that the underlying processes for the formation of 
penumbral filaments are the same as those for umbral dots. From 
this point of view, PUDs with dark-cored tails result from mag- 
netoconvective energy transport in the form of hot rising plumes 
inside the inclined magnetic field at the periphery of the um- 
bra. They represent a transition from CUDs to PGs and bright 
penumbral filaments. 

We conclude that the numerical MHD s imulations of mag- 
netoco nve ction in the sunspot umbra by ISchiissler & Voglerl 
(2006) and lRempel et alj d2009t) correctly predict and reproduce 
the observed substructures - dark lanes in CUDs and tails at- 
tached to PUDs. 
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